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Abstract 

Schroedinger equation on a Hilbert space 7i, represents a linear 
Hamiltonian dynamical system on the space of quantum pure states, 
the projective Hilbert space PTC. Separable states of a bipartite quan- 
tum system form a special submanifold of PTC. We analyze the Hamil- 
tonian dynamics that corresponds to the quantum system constrained 
on the manifold of separable states, using as an important example 
the system of two interacting qubits. The constraints introduce non- 
linearities which render the dynamics nontrivial. We show that the 
qualitative properties of the constrained dynamics clearly manifest the 
symmetry of the qubits system. In particular, if the quantum Hamil- 
ton's operator has not enough symmetry, the constrained dynamics is 
nonintegrable, and displays the typical features of a Hamiltonian dy- 
namical system with mixed phase space. Possible physical realizations 
of the separability constraints are discussed. 
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1 Introduction 



Classical and quantum descriptions of a physical system that is considered 
as composed of interacting subsystems have radically different features. The 
typical feature of quantum dynamics is the creation of specifically quantum 
correlations, the entanglement, among the subsystems. On the other hand, 
the typical property of classical description is the occurrence of chaotic or- 
bits and fractality of the phase space portrait, which can be considered as 
typically classical type of correlations between the subsystems. The type 
of correlations introduced by the dynamical entanglement does not occur in 
the classical description, and likewise, the type of correlations introduced 
by the chaotic orbits with fractal structures does not occur in the quantum 
description. This intriguing complementarity of the two descriptions repre- 
sents a problem that is expected to be solved by a detailed formulation of 
the correspondence principle. 

Comparison of typical features of classical and quantum mechanics is fa- 
cilitated if the same mathematical framework is used in both theories. It is 
well known, since the work of Kibble [lj,|2j,|3j, that the quantum evolution, 
determined by the linear Schroedinger equation, can be represented using the 
typical language of classical mechanics, that is as a Hamiltonian dynamical 
system on an appropriate phase space, given by the Hilbert space geometry 
of the quantum system. This line of research was later developed into the 
full geometric Hamiltonian representation of quantum mechanics. [1]-[12]. 
Such geometric formulation of quantum mechanics has recently inspired nat- 
ural definitions of measures of the entanglement [13], and has been used to 
model the spontaneous collapse of the state vector [S|,[I5], and dynamics of 
decoherence [IB] . 

It is our goal to use the geometric Hamiltonian formulation of quantum 
mechanics to study the relation between the dynamical entanglement and 
typical qualitative properties of Hamiltonian dynamics. Motivated by the 
fact that the Schroedinger equation can always be considered as a Hamil- 
tonian dynamical system, and that for Hamiltonian systems the definitions 
and properties of the dynamical chaos are well understood, we shall seek 
for a formal condition that when imposed on the Hamiltonian system repre- 
senting the Schroedinger equation of the compound quantum system renders 
the Hamiltonian dynamics nonintegrable and chaotic. It is well known that 
the linear Schroedinger equation of quantum mechanics represents always 
an integrable Hamiltonian dynamical system, irrespective of the dynamical 
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symmetries of the system. This is in sharp contrast with the Hamiltonian 
formulation of classical systems, where enough symmetry implies integrabil- 
ity and the lack of it implies the chaotic dynamics. Linearity of the quantum 
Hamiltonian dynamics, and the consequent integrability, is introduced in the 
Hamiltonian formulation by a very large dimensionality of the phase space of 
the quantum system. This high dimensionality can be considered as a conse- 
quence of two reasons. For a single quantum system, say a one dimensional 
particle in a potential, linear evolution and with it the principle of state 
superposition require infinite dimensional phase space of the Hamiltonian 
formulation. If the classical mechanical model is linear, say the harmonic os- 
cillator, the quantum Hamiltonian dynamics can be exactly describe on the 
reduced finite-dimensional phase space, the real plane in the case of the har- 
monic oscillator. The other related reason that increases the dimensionality 
of the quantum phase space compared to the classical model is in the way 
the state space of the compound systems are formed out of the components 
state spaces in the two theories. In order to represents the entangled states as 
points of the quantum phase space the dimensionality of the quantum phase 
space is much larger than just the sum of the dimensions of the components 
phase spaces. The points in the Cartesian product of the components phase 
spaces represent the separable quantum states and form a subset of the full 
quantum phase space. Needles to say, although the separable states are the 
most classical-like states of the compound system, they still are quantum 
states with nonclassical properties like nonzero dispersion of some subsys- 
tem's variables. Our main result will be that when the quantum dynamics, 
represented as a Hamiltonian system, is constrained on the manifold of sepa- 
rable quantum states the relation between the symmetry and the qualitative 
properties of the dynamics such as integrability or chaotic motion is reestab- 
lished. Thus, suppression of dynamical entanglement is enough to enable 
manifestations of the qualitative differences in dynamics of quantum systems 
and the relation between integrability and symmetry, traditionally related 
with classical mechanical models. 

In order to study the relation between the dynamical entanglement, sep- 
arability and the properties of Hamiltonian formulation of the quantum dy- 
namics we shall use, in this paper, the simplest quantum system that displays 
the dynamical entanglement, that is a system of two interacting qubits: 

H = luo- 1 + ucr 2 + n x OxOx + ii y Oy(jy + /i z <r z V 2 2 , (1) 

where cr x ,y,z l are the three Pauli matrices of the i-th qubit, and satisfy the 
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usual SU(2) commutation relations. In particular we shall compare the dy- 
namics of the system (1) in the case fi z ^ 0, fi x = fi y = with the case 
when \i x ^ 0, fi z — fj, y — 0. The former case is symmetric with respect to 
SO(2) rotations around z-axis and the later lacks this symmetry. Besides 
its simplicity, the systems of the form (1) are of considerable current interest 
because the hamiltonian of the universal quantum processor is of this form 

mm- 

Various lines of research, during the last decade, improved the under- 
standing of the relation between dynamical entanglement and properties of 
the dynamics. Strong impetus to the study of all aspects of quantum entan- 
glement came from the theory of quantum computation [18]. Quantization 



of classical non-integrable systems, and various characteristic properties of 
resulting quantum systems, have been studied for a long time [19]. The de- 
pendence of the dynamical entanglement, between a quantum system and its 
environment, on the qualitative properties of the dynamics of the system was 
studied indirectly, within the theory of environmental decoherence [20]. The 
relation between the rates of dynamical entanglement and the qualitative 
properties of the dynamics in the semi-classical regime was initiated in the 
reference [21] and various aspects of this relations have been studied since 
|22]-[31j. The relation between the symmetry of the genuinely quantum sys- 
tem (1) and the degree of dynamical entanglement was studied in reference 
[52] . As we shall see, our present analyzes is related to the quoted works, 
but the relation between the dynamical entanglement and symmetry is here 
approached from a very different angle 

The structure of the paper is as follows. We shall first recapitulate the 
necessary background such as: the complex symplectic and Riemannian ge- 
ometry of CP"; Hamiltonian formulation on CP™ of the quantum dynamics; 
geometric formulation of the set of separable pure states and Hamiltonian for- 
mulation of the constrained dynamics. In parallel with the general reminder, 
the explicit formulas for the system of two interacting qubits will be given. 
These are then applied, in section 3, to the study the qualitative properties 
of the separability constrained dynamics for the qubits systems. The main 
results are summarized and discussed in section 4. There we also discuss a 
model of an open quantum system with dynamics that clearly differentiates 
between the symmetric and the nonsymmetric systems. 
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2 Geometry of the state space CP 



Hamiltonian formulation of quantum mechanics is based on the fact that the 
scalar product of vectors \if) > in the Hilbert space of a quantum system can 
be used to represent the linear Schroedinger equation of quantum mechanics 
in the form of Hamilton's equations. The canonical phase space structure of 
this equations is determined by the imaginary part of the scalar product, and 
the Hamilton' s function is given by the quantum expectation < ip\H\ip > of 
the quantum hamiltonian. 

However, due to phase invariance and arbitrary normalization the proper 
space of pure quantum states is not the Hilbert space used to formulate the 
Schroedinger equation, but the projective Hilbert space which is the manifold 
to be used in the Hamiltonian formulation of quantum mechanics. In general, 
the resulting Hamiltonian dynamical system is infinite-dimensional, but we 
shall need the general definitions only for the case of quantum system with 
finite-dimensional Hilbert space, like the finite collection of qubits, in which 
case the quantum phase space is also finite-dimensional. We shall first review 
the definition of the complex projective space CP n , and then briefly state 
the basic definitions and recapitulate the formulas which are needed for the 
Hamiltonian formulation of the quantum dynamics on the state space and 
its restriction on the separable state subset. The general reference for the 
mathematical aspects of complex differential geometry is [33]. All concepts 
and formulas will be illustrated using the system of two interacting qubits. 

Differential geometry of the state space CP n is discussed by viewing it 
as a real 2n dimensional manifold endowed with complex, Riemannian and 
symplectic structure. In the case of CP™ this three structures are compatible. 

2.1 Definition and intrinsic coordinates of CP n 

States of a collection of N = n + 1 qubits are represented using normalized 
vectors of the complex Hilbert space C^. Since all quantum mechanical 
predictions are given in terms of the Hermitian scalar product on C^, and 
this is invariant under multiplication by a constant (vector independent) 
phase factor, the states of the quantum system are actually represented by 
equivalence classes of vectors in C^. Two vectors ip± and tp2 are equivalent: 
ip 2 ~ ipi if there is a complex scalar a ^ such that ip2 = aipi- This set of 
equivalence classes defines the complex projective space: CP n := (C n+1 — 
0)/ ~. It is the state space of the system of N qubits. Global coordinates 
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(c l ,...c N ) of a vector in that represent an equivalence class [ip], that 
is an element of CP™, are called homogeneous coordinates on CP n . The 
complex projective space is topologically equivalent to S' 2n+1 /S' 1 , where the 
In + 1-dimensional sphere comes from normalization and the circle S l takes 
care of the unimportant overall phase factor. 

The projective space CP™ is locally homeomorphic with C n . Intrinsic 
coordinates on CP n are introduced as follows. A chart consists of equiv- 
alence classes of all vectors in (C ra+1 — 0) such that c M 7^ 0. In the chart 
the local ( so called inhomogeneous) coordinates ( u , v — 1, 2 . . . n are given 
by: 

C = C (!/</*- 1), C = C +1 (1/ > (2) 

where 

= v= l,2,.../i-l,/i+l,...n + l. (3) 

The coordinates C^( c ) an d C^'( c ) °f a point c which belongs to the domain 
where two charts and overlap are related by the following holomorphic 
transformation 

Q(c) = (c7^')c;(c) (4) 

As an illustration consider the system of two qubits. The Hilbert space 
is H = Hi<S>H2 = C 2 ® C 2 = C 4 . As a basis we can choose the set 
of separable vectors | Tt>> I tl>> I lt>> I ll> or an y other four orthogonal 
vectors. The coordinates of a vector in C 4 with respect to a basis are denoted 
(c 1 , c 2 , c 3 , c A ) . The corresponding projective space is CP 3 = S 7 /S 1 . At 
least two charts are needed to define the intrinsic coordinates over all CP 3 . 
Consider first all vectors with a nonzero component along |1 >= | j"f> that 
is c 1 7^ 0, i.e. all vectors except the vector | ||>. Then the numbers ^ 
are defined as = c 1 /^ = 1, £ 2 = c 2 /c 4 , £f = c'/c 1 , £ 4 = c^/c 1 and 
finally the three intrinsic coordinates (Ci,Ci>Ci) are gi ven by relabelling of 
£i : Ci 1 = £1, Ci = £i> Ci = £i- To coordinatize the vector |4 >= | ||> we 
need another chart. 

Quantum mean values of linear operators on C 4 are indeed reduced to 
functions on CP 3 . For example, consider the following Hamiltonian operator 

H = ujiy z ® 1 + cul <S> o z + ® &x (5) 

In the separable bases the normalized quantum expectation < ip\H\ip > / < 
ip\ip > is given by the following function of (c 1 , c 2 . . . , c 4 ) 

= 2c(c 1 c 1 - c 4 c 4 ) + Mc 2 c 3 + c 3 c 2 + cV + cV) 

c 1 ^ + c 2 c 2 + c 3 c 3 + c 4 c 4 ' 1 J 
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In the intrinsic coordinates C'SC^C 3 and their conjugates this expression is 
given by 

2c(i- C 3 C 3 )+MC 1 C 2 + C 2 C 1 + C 3 + C 3 ) 
1 + ec 1 + C 2 C 2 + C 3 C 3 

We shall also analyze the following Hamiltonian 



H = ucr z ® 1 + uj! ® a z + fia z <g> cr z , (8) 
whose normalized mean value is given by 

2uj(c 1 c 1 - c 4 c 4 ) + fijc'c 1 + c 4 c 4 - c 2 c 2 - C 3 C 3 ) 

c lcl + C 2c2 + C 3c3 + C 4c4 " ^ ^ 

The corresponding function on CP 3 is, in the intrinsic coordinates, given by 

C 3 C 3 ) + Mi + C 3 C 3 - - C 2 C 2 ) nm 

i + ec 1 + c 2 c 2 + c 3 c 3 ' 1 ) 

2.1.1 Submanifold of separable states 

Consider two quantum systems A and B with the corresponding Hilbert 
spaces Ha and Hb- Taken together, the systems A and B form another 
quantum system. The statistics of measurements that could be performed 
on this compound system requires that the Hilbert space of the compound 
system is given by the direct product Hab = Ha ®Hb- The space of pure 
states of the compound system is the projective Hilbert space PHab- in the 
case of finite dimensional state spaces PH n A +l = CP n and PH n A l+l = CP m the 
state space of the compound system is (7p( m + 1 )(™+ 1 )- 1 . Vectors in Hab of the 
form ipA <8> i>B where ipA/B £ Ha/b are called separable. The corresponding 
separable states form the (m + n)-dimensional submanifold CP m x CP n 
embedded in cp(™+i)(™+i)-i. 

In the case of two qubits the submanifold of the separable states CP 1 x 
CP 1 forms a quadric in the full state space CP 3 , given in terms of the 
homogeneous coordinates (c 1 , c 2 , c 3 , c A ) of CP 3 by the following formula 

cV = cV. (11) 

In terms of the intrinsic coordinates C 1 ? C 2 ? C 3 ? i n the chart with c 1 ^ 0, i.e. 
= 1, the equation (11) is 

C X C 2 = C 3 . (12) 
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2.2 Complex structure on CP n 

Consider a complex manifold M. with complex dimension dim c M. — n (in 
particular CP™ ). We can look at M. as a real manifold with (\xm R M = 2n. 
The real coordinates (x 1 , . . . x 2n ) are related to the holomorphic • • • C n ) 
and anti-holomorphic (C 1 , . • • C n ) coordinates via the following formulas: 

(x" + ix" +n )/V2 = C v = 1,2,... n, 

{x v -ix u+n )/V2 = C, v = 1,2,... n, (13) 

and 



f = x v = (C + 0/V2, I / = l,2,...n, 
p^ = ^ + " = (C"-r)/A i/ = l,2,...n. (14) 



The tangent space T X A^ is spanned by 2n vectors: 

9 d d d 

or by the basis 

d d d d 

An almost complex structure on a real 2n-dimensional manifold is given 
by a (1, 1) tensor J satisfying J 2 = 1, i.e. J" J£ = — 5 a . Locally, the almost 
complex structure J is given in the real coordinates by the following matrix 



-1 

1 



(17) 



where 1 is n-dimensional unit matrix. If the real In manifold is actually a 
complex manifold, like in our case, the almost complex structure is defined 
globally and is called the complex structure. 

2.3 Riemannian structure on CP n 

Hermitian scalar product induces a complex Euclidean metric on C^. The 
metric induced on CP™ is the Fubini-Study metric, and is given, in (£, () 
coordinates, using annxn matrix with following entries 

VivO ^ = U...n, (18) 
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where CC = EpC*e- 

The Fubini-Study metric in (C,C) coordinates is then given by 2n x 2n 
matrix 

In the real coordinates the Fubini-Study metric is given by the standard 
transformation formulas 



Gij(Q,P) = gfc,i(C(g»P),C(g,P)) g^* » 

where we used Z = . . . ( n and X = (q 1 . . .p n ). 

In the example of two qubits the Fubini-Study metric on CP 3 is 



(20) 







(l+CO-^C 1 
(i+CC) 2 



(i+CC) 2 

V (i+cc) 2 







-c 2 c -1 

(i+CC) 2 
(i+CC)-C 2 C 2 

(I+CC) 2 
-c 2 c 3 



2G 





(i+CC) 2 

-c 3 c 2 

(i+CC) 2 
(i+CC)-C 3 C 3 
(i+CC) 2 



(l+CO-C 1 ^ 


_ C l C -2 


_ C l f -3 


(i+CC) 2 

_ C 2 C "1 


(i+CC) 2 
(i+C<)-C 2 C 2 


(i+CC) 2 

-c 2 c 3 


(i+CC) 2 

-c 3 c _1 


(i+CC) 2 

-c 3 c 2 


(i+CC) 2 
(i+CC)-C 3 C 


(i+CC) 2 


(i+CC) 2 


(i+CC) 2 






























(i+CC) 2 (i+CC) 2 w w w , / 

(21) 

Transformation to the real coordinates, by application of the formula (20), 
gives 



V 



a 

p 1 p 2 +q 1 q 2 



pV-pV 

a 

p 1 q 3 —p 3 q 1 



p 2 p 3 +q 2 q 3 
a 



p 2 q 3 — q 2 p 3 



a 

p 2 p 3 +q 2 q 3 







p a q 1 —p 1 q 3 
a 

p 3 q 2 -p 2 q 3 



a 





a 
b 

p 1 p 2 +q 1 q 2 
a 

p^p 3 — q 1 q 3 



19 9 1 

p q —p q 
a 



p 3 q 2 —p 2 q 3 
a 

p 1 p 2 +q 1 q 2 



a 

p 2 p 3 +q 2 q 3 



p 1 q 3 —p 3 q 1 \ 
a 

p 2 q 3 —p 3 q 2 



_ p 1 p 3 q 1 q 3 
a 

p 2 p 3 +q 2 q 3 



where 

a = (p 1 ) 2 +(p 2 ) 2 +( P 3 ) 2 +(g 1 ) 2 +(g 2 ) 2 +(g 3 ) 2 +2, b 
Obviously, Q is positive definite and symmetric. 



(22) 

(p 1 ) 2 +(p 3 ) 2 +(g 1 ) 2 +(g 3 ) 2 +2- 
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2.4 Symplectic structure on CP n 

The Hermitian scalar product on is also used to define the symplectic 
structure on C N and this induces the symplectic structure on CP™. The 
symplectic structure is the closed nondegenerate two form Q on CP™, which 
is, in (£, () coordinates given by 

" = «0(C, CWC" A C" (23) 

where is the Fubini-Study metric (18). In real coordinates, the symplectic 
structure is given by fl(q,p) = JG(q,p) where G(q,p) is given by (20) and J 
by (17). 

The symplectic form on the two qubits state space is in the real bases 
given by the product of matrices (17) and (22). The results is 







n = 










pV-pV 


— p 2 q 1 +p 1 q 2 


—p i q 1 +p 1 q 3 


b 


p 1 p 2 +q 1 q 2 


p 1 p 3 +q 1 q s \ 


a 2 



pV-pV 


a 2 

p' 2 q i —p i q 2 


a 2 

p 2 p 1 +g 1 g 2 


a 2 
b 


a 2 

p 2 p 3 +q 2 q 3 


a 2 

p i q 1 -p 1 q i 


a 2 



P 1 P 3 +g 1 <? 3 


a 2 

p 1 p 3 + (J 1 Q 3 


a 2 

p 2 p :i +q 2 q :i 


a 2 
b 


a 2 
b 


a 2 

p 2 p 1 +q 1 q 2 


a 2 



p 2 q 1 —p 1 q 2 


a 2 

p 1 q 2 —p 2 q 1 


a 2 

p 1 q 3 —p 3 q 1 


a 2 


a 2 
b 


a 2 


a 2 



p :i q 2 '—p 2 q 3 


a 2 

p 2 q 3 —p 3 q 2 


a 2 

p 1 p 3 +q 1 q 3 


a 2 

p 2 p :i +q 2 q 3 


a 2 
b 


a 2 

p 3 ^ 1 — p 1 ^ 3 


a 2 

o J 


a 2 


a 2 


a 2 


a 2 


a 2 



(24) 



3 Quantum Hamiltonian dynamical system 
on CP n 

The Schroedinger equation on C N is in some basis >, % = 1,2. . . N} 
given by: 

i—^<^\H\^>d>. (25) 

In the real coordinates this equation assumes the form of a Hamiltonian 
dynamical system on R 2N with a global gauge symmetry corresponding to 
the invariance \ip >— > exp(ix)\ip >. Reduction with respect to this symmetry 
results in the Hamiltonian system on CP™, considered as a real manifold with 
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the symplectic structure given by (23). The Hamilton equation on CP n , that 
are equivalent to the Schroedinger equation (25), are 

a J. 

^ = 2n l > k V k H(x), (26) 

where Q l ' k is the inverse of the symplectic form, and H(x) is given by the 
normalized quantum expectation of the Hamilton's operator < ip\H\ip > 
I < > expressed in terms of the real coordinates (14). For example, 
the hamiltonian (7) is given in terms of the real coordinates q l = x l , p % = 
x l+n , i — 1, . . . n by 

H = -[2- (p 3 ) 2 - (q 3 ) 2 ] + %V + qV + V2q 3 ). (27) 
a a 

and the symmetric hamiltonian (9) is given by 
H = -^-(p 3 ) 2 -^) 2 ]-^^ 1 )^^)^^)^^) 2 -^) 2 -^) 2 -^ (28) 

(X (X 

The Hamilton's equations (26) with the hamiltonian (27) and the sym- 
plectic form (24) assume the following form 

q 1 = -2ujp l + fip 2 - /i(pV +p 1 q 3 )/V2 

q 2 = -2up 2 + fip 2 - fi(p 3 q 2 +p 2 q 3 )/V2 

q 3 = -Aujp 3 - V2/ip 3 q 3 

p 1 = 2cjg 1 -/ig 2 + /i(gV-pV)/ v / 2 

p 2 = 2uq 2 - fiq 1 + fi(q 3 q 2 -p 2 p 3 )/V2 

p 3 = 4uq 3 + ^(q 3 ) 2 -(p 3 ) 2 -2)/V2. (29) 

The equations of motion with the symmetric hamiltonian (28) on CP 3 
are quite simple 

q 1 = -2(uj + fj^p 1 

q 2 = -2(u + fx)p 2 

q 2 = -Aup 3 

p 1 = 2{u + ^)q l 

p 2 = 2(to + fi)q 1 

p 3 = Aujq 3 . (30) 
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3.1 Quantum Hamiltonian system with imposed sep- 
arability constraints 

Dynamics of a constrained Hamiltonian system is usually described by the 
method of Lagrange multipliers [34] . |35] . Consider a Hamiltonian system 
given by a symplectic manifold M. with the symplectic form Q and the 
Hamilton's function H on M. . Suppose that besides the forces described 
by H the dynamics of the system is affected also by forces whose sole effect 
is to constrain the motion on a submanifold TV G M. determined by a set 
functional relations 

fi(q,p) = ---fk(q,p) = (31) 

The method of Lagrange multiplies assumes that the dynamics on M is de- 
termined by the following set of differential equations 

k 

X = fi(VX, VH'), H' = H + J2, *jfj (32) 

j 

which should be solved together with the equations of the constraints (31). 
The Lagrange multipliers Xj are functions of (p, q) that are to be determined 
from the following, so called compatibility, conditions. 

f^QiVfaVH') (33) 

on M. The equations (33) uniquely determine the functions Ai(p, <?),... \k(p, l) 
if and only if the matrix of Poison brackets {fi, fj} = Q(Vfi, V/j) is nonsin- 
gular. If this is the case then all constraints (31) are called primary, and M 
is symplectic manifold with the symplectic structure determined by the so 
called Dirac-Poison brackets 

{F u F 2 }' = {F u F 2 } + F 1 }{f l , f 3 Y l {f v F 2 } (34) 

As we shall see, this is the case in the examples of pairs of interacting qubits 
constrained on the manifold of separable states that we shall analyze. On the 
other hand, if some of the compatibility equations do not contain multipliers, 
than for that constrain fj = {fj,H} = 0, which represents an additional 
constraint. These are called secondary constraints, and they must be added 
to the system of original constraints (31). If this enlarged set of constraints 
is functionally independent one can repeat the procedure. At the end one 
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either obtains a contradiction, in which case the original problem has no 
solution, or one obtains appropriate multipliers A& such that the system (33) 
is compatible. In the later case the solution for might not be unique in 
which case the orbits of (32) and (31) are not uniquely determined by the 
initial conditions. 

Let us apply the formalism of Lagrange multipliers on the system of 
two interacting qubits additionally constrained to remain on the manifold 
of separable pure state. The real and imaginary parts of (12) give the two 
constraints in terms of real coordinates {q l ,q 2 ,q 3 ,p l ,p 2 ,p 3 ) 

fi = pY - q X q 2 + V2q 3 , f 2 = V2q 3 - p 2 q l - p l q 2 (35) 

The compatibility conditions (33) assume the following form 

A = ft(V/i,Vif) + A 2 ft(V/i,V/ 2 ) = 0, 

h = ft(V/ 2 ,Vif) + Aift(V/ 2 ,V/i) = 0. (36) 

where VL is the symplectic form (24) and ft(V/i, VH) = tt a > b V a fiVbH. 
The matrix of Poisson brackets {fi, fj} on Af is 

[2 + (p 1 ) 2 + {q l f}[2 + (p 2 ) 2 + (q 2 ) 2 ]/8 

-[2 + (p 1 ) 2 + (q^m + (p 2 ) 2 + (q 2 ) 2 ]/8 

(37) 

and is nonsingular. Thus the compatibility conditions can be solved for the 
Lagrange multipliers Xi(q,p), A 2 (g,p), 

Ai = 4 VpW + [(g 1 ) 2 - 2] [2 + (p 2 ) 2 - (g 2 ) 2 } + (p 1 ) 2 ^) 2 - (p 2 ) 2 - 2] 



A 2 = S/j, 



[2 + (p 1 ) 2 + (q 1 ) 2 ) 2 (2 + (p 2 ) 2 + (q 2 ) 2 ] 2 
(p 1 ) 2 p 2 q 2 - p 2 q 2 [(q 1 ) 2 - 2] + p l q l \2 + (p 2 ) 2 - (q 2 ) 2 ] 
[2 + (p 1 ) 2 + {q 1 ) 2 ) 2 (2 + (p 2 ) 2 + (q 2 ) 2 ] 2 

(38) 



Finally, the dynamics of the constrained system is described by the equa- 
tions (32) and (31) with Xi(q,p), A 2 (q , ,p) and fi(q,p), f2(q,p) given by (38) 
and (35). For the Hamiltonian (27) the resulting equations of motion for 
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q , q ,p ,p are 



q 



1 Ajip 1 q 1 q 2 + 2UJP 1 [2 + (p 2 ) 2 + (q 2 ) 2 } 



2 + (p 2 ) 2 + (q 2 ) 2 
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. 2 4jjpW - 2up 2 [2 + (p 1 ) 2 + (g 1 ) 2 ] 

9 ~ 2 + (p 1 ) 2 + (g 1 ) 2 

.! 2^g 2 [(g 1 ) 2 - (p 1 ) 2 - 2] + 2^2 + (p 2 ) 2 + (g 2 ) 2 ] 
P 2 + (p 2 ) 2 + (g 2 ) 2 

2/ig 1 [(g 2 ) 2 - (p 2 ) 2 - 2] + 2cg 2 [2 + (p 1 ) 2 + (g 1 ) 2 ] 
P " 2 + (p 1 ) 2 + (g 1 ) 2 • 1 ' 

The same procedure for the symmetric hamiltonian (28) results with the 
following equations of motion 



p 1 



2/ip 1 [(P 2 ) 2 + (g 2 ) 2 - 2)] - 2oy[2 + (p 2 ) 2 + (g 2 ) 2 ] 

2 + (p 2 ) 2 + (g 2 ) 2 
2/xp 2 [(p 1 ) 2 + (g 1 ) 2 - 2] - 2a;p 2 [(2 + (p 1 ) 2 + (g 1 ) 2 ] 

2+ (pl)2 + (gl)2 

-2/ig 1 [(g 2 ) 2 + (p 2 ) 2 - 2] + 2<V[2 + (p 2 ) 2 + (g 2 ) 2 ] 

2 + (p 2 ) 2 + (g 2 ) 2 
-2/ig 2 [(g 1 ) 2 + (p 1 ) 2 - 2] + 2cug 2 [2 + (p 1 ) 2 + (g 1 ) 21 



p2 ^ 2 + (p 1 ) 2 + (g 1 ) 2 ' ' ' (40) 

There are also the equations expressing g 3 and p 3 in terms of g 1 , g 2 '^p 1 ,p 2 , 
but the solutions of these are already given by the constraints. 



3.2 Qualitative properties of the constrained dynam- 
ics of two interacting qubits 

In this section we present the results of numerical analyzes of the qualita- 
tive properties of the dynamics generated by the constrained equations (40) 
and (39), corresponding to the quantum Hamiltonians (28) with the 5*0(2) 
symmetry and (27) without such symmetry. 

It is well known that any quantum system is integrable when considered 
as the Hamiltonian dynamical system on the symplectic space 7i, and that 
the reduction on the symplectic manifold P7i preserves this property. This is 
simply a consequence of the form of the quantum Hamiloton's function, which 
is always defined as the mean value of the Hamiltonian operator. Contrary 
to the case of classical Hamiltonian systems, the symmetry of the physical 
system has no relevance for the property of integrability in the Hamiltonian 
formulation of the Schroedinger equation. We illustrate this fact, in figures 
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Figure 1: Projections on (qi,pi) plane of a typical orbit for the hamiltonian 
systems (28) (a) and (27) (b) on CP 3 and on the submanifold of separable 
states (c) for (40) and d for (39). The values of the parameters are uj — 1 
and /i — 1.7 

la,b, by projections on (q 1 ,^ 1 ) plane of a typical orbit for the symmetric and 
nonsymmetric hamiltonians of the pair of qubits. The motion on CP 3 in 
the symmetric case has further degeneracy compared with the nonsymmetric 
case, but both cases generate integrable, regular Hamiltonian dynamics. 

On the other hand, the qualitative properties of the dynamics constrained 
by the separability conditions, are quite different. Typical orbits in the sym- 
metric and nonsymmetric cases are illustrated in figure lc,d. Symmetric 
dynamics constrained by separability is still regular, while the nonsymmetric 
Hamiltonian generates the constrained dynamics with typical chaotic orbits. 
This is further illustrated in figures 2, where we show Poincare surfaces of 
section, defined by q 2 = 0,p 2 > and H(p l ) q 1 , p 2 , q 2 ) = h for different values 
of the coupling //. Obviously, the constrained system displays the transition 
from predominantly regular to predominantly chaotic dynamics, with all the 
intricate structure of the phase portrait, characteristic for typical Hamilto- 
nian dynamical systems. Thus, we can conclude that the quantum system 
constrained on the manifold of separable state behaves as typical classical 
Hamiltonian systems. If there is enough symmetry, i.e. enough integrals of 
motion, the constrained dynamics is integrable, otherwise the constrained 
quantum dynamics is that of typical chaotic Hamiltonian system. 
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Figure 2: Poincare sections for the separability constrained non-symmetric 
quantum dynamics (39). The parameters are to — 1, h — 1.5 and (a) /i = 1.3, 
(b) n = 1.7 



4 Summary and discussion 

We have studied Hamiltonian formulation of quantum dynamics of two in- 
teracting qubits. Hamiltonian dynamical system on the state space CP 3 
as the phase space, is integrable irrespective of the different symmetries of 
the quantum system. We have then studied the dynamics of the quantum 
Hamiltonian system constrained on the manifold of separable states. The 
main result of this analyzes, and of the paper, is that the quantum Hamilto- 
nian system without symmetry generates nonintegrable chaotic dynamics on 
the set of separable states, while the constrained symmetric dynamics gives 
an integrable system. It is important to bare on mind that neither the system 
nor the separable states that lie on an orbit of the constrained system have 
an underlining classical mechanical model. Thus, forcing a non-degenerate 
quantum system to remain on the manifold of separable states is enough to 
generate a dynamical system with typical properties of Hamiltonian chaos. 

Our analyzes of the separability constrained quantum dynamics has been 
rather formal. In order to inquire into possible interpretation of our results we 
need a model of a physical realization of the separability constraints. To this 
end we consider an open quantum system of two interacting qubits, whose 
dynamics satisfies the Markov assumption [36], and we choose a Hermitian 
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Lindblad operator of the following form 

L = Zncr+eri <g> a 2 + a 2 _ + li 2 a\a l _ <g> a 2 _o\ + l 2 \o\a\ ® a 2 + a 2 _ + X % iO X _o\ <8> o 2 _o\ 

2 

= kj\i Xil 1 ® \i ><j\ 2 (41) 

where |1 >= | f> and |2 >= | |>. 

The dynamics of a pure state of the open system under the action of a 
Hamiltonian H and the Linblad 7L is described by the following stochastic 
nonlinear Schroedinger equation [36].[37] 

|# > = -iH\if> >dt+ -j{L- < ip\L\ip >) 2 \ip > dt 

+ 7(L- < >)|^ > dW, (42) 

where (iW 7 is the increment of complex Wiener c-number process W(t). 

The equation (42) represent a diffusion process on a complex Hilbert 
space, and is central in the "Quantum State Diffusion" (QSD) theory of 
open quantum systems [37j. It has been used to study the systems of in- 
teracting qubits in various environments for example in [32], [38], and the 
effect of the Linblad operator (41) on the entanglement between two qubits 
was considered in [16]. The influence of the non-Hamiltonian terms of drift 
(proportional to 7 2 ) and the diffusion (proportional to 7)), with the Linblad 
operator of the form (41), is to drive an entangled state towards one of the 
separable states with the corresponding probability. This process occurs on 
the time scale proportional to 7 _1 . So, for large 7 there occurs an almost 
instantaneous collapse of an entangled state into a separable one. We believe 
that with a proper choice of the parameters lij the long term dynamics of 
a pure state described by (42) can have the same qualitative properties as 
the separability constrained quantum dynamics. In particular, the difference 
between the qualitative properties of symmetric and nonsymmetric systems, 
reflected in the constrained Hamiltonian system, should also manifest in the 
dynamics of (42) for a proper choice of kj. This expectations are supported 
by figures 3, which illustrate the dynamics of (< a], >, < a 2 >) for the Hamil- 
tonian operators (5) and (8) as calculated using the constrained Hamiltonian 
equations (39) and (40) (figures 3b and 3a ), or the QSD equation (42) (fig- 
ures 3d and 3c) for a particular choice of kj and large 7 = 5. Of course, the 
choice of optimal values for should be according to some criterion, which 
is the problem we are currently investigating. 
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Figure 3: Figures illustrate the dynamics of (< a x >,< a y >) for the 
constrained Hamiltonian systems (40) (a) and (39) (b) and for the stochastic 
Schroedinger equation (42) with the Linblad (41) and the hamiltonians (8) (c) 
and (5) (d). the parameters are u — 1, // = 1.7, 7 = 5 and = 0.21, 11, 2 = 
0.21, Z 2 ,i = 0.215,^,2 = 0.205. 

The pair of coupled qubits, analyzed in this paper, is the simplest quan- 
tum system exhibiting dynamical entanglement. We intend to investigate 
the effects of suppression of the dynamical entanglement in systems with 
spacial degrees of freedom, obtained by quantization of classically chaotic 
systems, for example a pair of coupled nonlinear oscillators. In this case, 
the Hamiltonian formulation of the quantum dynamics requires an infinite- 
dimensional phase space, and the analyzes of the separability constrained 
dynamics is more complicated. However, it wold be interesting to compare 
the dynamics obtained by separability constraints with that of some more 
standard semi-classical approximation. 
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FIGURE CAPTIONS 

Figure 1 Projections on (qi,pi) plane of a typical orbit for the hamil- 
tonian systems (28) (a) and (27) (b) on CP 3 and on the submanifold of 
separable states (c) for (40) and d for (39). The values of the parameters are 
uj = 1 and n = 1.7 

Figure 2 Poincare sections for the separability constrained non-symmetric 
quantum dynamics (39). The parameters are ou — 1, h — 1.5 and (a)// = 1.1, 
(b) \x = 1.3, (c) fi = 1.5 and (d) fi = 1.7 

Figure 3 Figures illustrate the dynamics of (< a x >,< cr y >) for the 
constrained Hamiltonian systems (40) (a) and (39) (b) and for the stochastic 
Schroedinger equation (42) with the Linblad (41) and the hamiltonians (8) (c) 
and (5) (d). the parameters are cu — 1, (j, — 1.7, 7 = 5 and = 0.21, 11, 2 = 
0.21, l 2 ,i = 0.215,^,2 = 0.205. 
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